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teracting bosons in the Standard Model and in model theories is discussed. In 
particular, such processes are associated with the baryon and lepton number vio- 
lation in the Standard Model. The most interesting domain where the multiplicity 
of bosons is larger than the inverse of small coupling constant is beyond the scope 
of perturbation theory and requires a non-perturbative analysis. 
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1 Introduction 



The vastness of the field of non-perturbative methods in high-energy physics inevitably 
compels me to focus on a specific topic among those which attract a considerable in- 
terest and where a non-trivial development is likely in the near future. One such topic, 
which also has dominated the parallel session on non-perturbative methods, is related to 
multiboson phenomena in the electroweak physics and, more generally, in models with 
weak coupling. These phenomena are interesting because of two basic reasons. One is 
that it is with multiboson processes is associated the violation of the sum of the baryon 
(B) and the lepton (L) numbers in the Standard Model]]]]. Therefore such processes 
determine the evolution of (B+L) at high temperature in the early universe 0]. Also as 
initially envisioned in early works || || and indicated by specific calculations]^, |(| the 
processes with (B+L) violation and production of many electroweak bosons might be in 
principle observable in high energy collisions in the multi-TeV energy range. The other, 
theoretical, reason is related to the old-standing problem of the factorial divergence of 
perturbation theory series, which dates back to the work of Dyson [0. This problem looks 
to be a matter of a purely theoretical concern as long as the quantities under discussion 
are such that they appear at low orders, like the anomalous magnetic moment. For such 
quantities the inability in principle to find the exact result by the perturbative expan- 
sion, though disappointing, does not prevent from calculating in few first orders with an 
accuracy required by practical measurements or greater. However the problem of the 
inherent divergence of perturbative series becomes quite acute as soon as one considers 
processes at energies such that a large number of interacting particles can be produced, 
i.e. the processes which occur starting only from a high order of the perturbation theory, 
where the expansion becomes unreliable. 

At present there is a general understanding that multiparticle electroweak processes 
with many bosons both in the initial and the final state {many — > many scattering), 
including those with (B+L) violation, are not suppressed at high temperature and thus 
they indeed determine the (B+L) history of the universe. On the other hand the un- 
derstanding of the processes, in which many bosons are produced by two or few initial 
particles with high energy {few — > many scattering), is far from complete and there 
are arguments both pro and con the idea that at sufficiently high multiplicity of final 
particles such processes can have an observable cross section. 

The contribution of the many — > many scattering at high temperature is described 
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within the WKB technique by expansion around non-trivial classical solutions of the 
field equations: sphalerons|8], [| and, more generally, periodic instantons[|TO[|. It is also 
believed that the few —>■ many scattering can be fully described by applying a WKB 
technique using special classical configurations of the field. However, this is still a 
conjecture, and a specific method of a full WKB analysis of the latter scattering has yet 
to be developed. 



2 (B+L) violating electroweak processes 

As is known|l[], the electroweak interaction in the Standard Model does not conserve 
the sum of the baryon and the lepton numbers as a result of the triangle anomaly. The 
amount of the (B+L) violation in a process is determined by the change of the winding 
number N C s of the electroweak gauge fields: 

A(B + L) = 6AN CS . (1) 

However, changing the winding number by one or several units requires presence (at 
least in the intermediate state) of the W and Z field configurations with energy of order 
my/ /aw- This is usually illustrated by the sketch of the dependence of minimal energy 
of the field with a given Nqs as a function of Ncs shown in figure 1, and the field 
configuration, corresponding to the top of the barrier, is the so-called sphaleron|| with 
energy E$ p about 10 TeV [[|. 

If the energy E available in a process is much less than E$ p then the only way in which 
(B+L) can be violated is due to quantum tunneling, which at E = is described by 
the instanton[|TTJ solution to the Euclidean field equations, whose action is Si = 2n/aw- 



The amplitude of such process then contains the WKB tunneling factor exp(— 27r/«vy) ~ 
1CT 80 , which thus makes the process unobservable by any practical measure. 

The plot in the figure 1 however invites the suggestion that once the available energy 
is close to or larger than E$ p the suppression of (B+L) violating processes should weaken 
or disappear altogether. The two relevant situations where large energy is available in 
individual processes are high temperatures and high-energy particle collisions. 



2.1 (B+L) violation at high temperature 

As first realized by Kuzmin, Rubakov and Shaposhnikov@] the rate with which the 
system traverses the sphaleron barrier in thermal equilibrium at a temperature T < E$ p 
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Figure 1: The dependence of minimal energy of the electroweak gauge field on its winding 
number. 



is determined by the Boltzmann factor exp(— E Sp (T)/T) and may become unsuppressed 
at temperatures larger than E$ p . The dependence of E$ p on the temperature arises 
through the temperature dependence of the vacuum expectation value of the Higgs 
field, which sets the electroweak energy scale. In particular the v.e.v. vanishes at the 
phase transition temperature T c , above which the electroweak symmetry is restored. 
Thus at T > T c the sphaleron barrier is absent and the (B+L) violating processes may 
go without an exponential suppression. As a result [10, 12 - 14] the rate of change of 
the (B+L) density is given by 



A(B+L) 



dexpt-EspCH/T) if T<T C 
C 2 aLT 



if T > T r 



(2) 



where C\ and C*2 are constants. 

In particular, the prefactor C\ for the 'low' temperature rate is determined by fluc- 
tuations of the fields near the sphaleron configuration. The initial calculations [15 - 17], 
which considered only the contribution of bosonic fluctuations, were most recently rean- 
alyzed and extended [|l^] to include also the fermionic determinant. The contribution of 
the heavy top quark is found to significantly suppress the prefactor C\, which enables to 
somewhat relax the upper bound on the mass of Higgs boson in the minimal Standard 



Model, following from the requirement [H], [14], |19| that the (B+L) violating processes in 
the early universe immediately after the electroweak phase transition (i.e. just below T c ) 
do not wash out completely the baryon asymmetry, independently of the mechanism by 
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which it was created before or during the phase transition. Using their result for G\ and 
m t = 174 GeV, Diakonov et.al. [|18| find this upper bound to be m H < 66 GeV, which 



is only slightly higher than the lower bound mn > 58.4GeV [£(J from a direct search 
at LEP. Therefore an improvement in the experimental search can either find the Higgs 
boson or close the gap of compatibility of the minimal model with the observed baryon 
asymmetry of the universe. 

The shape of the sphaleron barrier in the presence of heavy top and the evolution 
of the energy levels of a heavy fermion were considered in detail in the contributed 



papers [21] and p2| respectively. Multi-sphaleron configurations are considered in [23 



and electroweak strings, viewed as "stretched sphalerons" in [24]. However the role of 



the latter configurations in thermal equilibrium is yet to be clarified. 
2.2 (B+L) violation in high-energy collisions 

The sphaleron energy scale E$ p is within (hopefully) reachable energies at prospective 
colliders. Therefore a most intriguing question arises as whether the exponential sup- 
pression of the (B+L) violating processes vanishes at an energy of order E$ p in collision 
of two leptons or quarks. The difference between the high-temperature (B+L) violation 
and the processes induced by just two or few energetic particles is that in the former case 
the dominant contribution to the rate comes [HJ from processes in the thermal bath, in 
which many soft particles with total energy E > E$ p scatter into a final state of also 
soft particles with different (B+L), while in the latter case a coupling between the hard 
initial particles and soft modes of the field with a non-trivial topology is required. 

Following the conjecture [|13|] that an enhancement of the cross section of (B+L) 
violating scattering may be associated with multiparticle final states, Ringwald|| and 
Espinosa|| pursued a calculation of a generic instanton-induced process of the type 

J + J^10f + n w W + n H H , (3) 

where nw (^h) is the multiplicity of produced gauge (Higgs) bosons and / stands for a 
quark or a lepton. (The presence in the instanton-induced scattering of twelve fermions: 
nine quarks and three leptons, one from each electroweak doublet is mandated by the 
anomaly condition in eq. ([!]), i.e. by the number of fermionic zero modes of an instanton.) 
The amplitude of the scattering (H) was found to factorially depend on the multiplicity 
of bosons: A ~ nw^n-H^ exp(— Si), which lead to the argument ||25|| that the factorial 



enhancement may beat the exponential suppression at riw,H > 0(1 /caw) i-e. at energy 
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larger than 0(Esp). In fact the growth with energy of the total cross section for (B+L) 
violating processes, observed in the early calculations [[5], ^, |25|, suggested p5| that this 
cross section may become strong: reach its unitarity limit at energies in the multi-TeV 
range. 



2.2.1 "Holy grail" function. 



By quite general scaling argument s [[26| , |27| the total cross-section of instanton-induced 
scattering should obey the scaling behavior 

.±L F (E 

Tnw/otw- The function F(e) is often termed as "holy grail" function. 



tot 
a A(_B+L) 



cxp 



(4) 



where E n ~ E\ 



s P 



At e = one has F(0) = 1, while the initial enhancement]^, |5| of the cross section due 
to opening multi-boson channels corresponds^, f29f to the first non-trivial term in the 
expansion in e: F(e) = 1 — §e 4//3 + • • ., where e = E/E with E = \fQnmw/oLw ~ 
18TeV. The expansion in fact goes in powers of e 2//3 , and by now two next terms are 
known [30 - 34]: 



F(e) 



e 8/3 ln £ + 



(5) 



o 16 32 \ m w , 

The latter two terms are determined by interaction between soft final particles. Starting 
from the term of order e 10 / 3 the "holy grail" function is also contributed by interactions 
between hard initial and soft final particles and by interaction between the initial hard 
particles [35 - 37]. 

Unfortunately, any finite number of terms in the expansion of F(e) does not allow 
to assert the behavior of the function at finite e ~ 0(1). Therefore it is not known yet, 
whether the function F(e): 

i) goes to zero at finite e (finite energy), 

ii) goes to zero as e — » oo, or 

in) is bounded from below by a positive value. 

Certainly, the most interesting phenomenologically is the first possibility, since then the 
cross section with (B+L) violation and multiboson production becomes observably large 
at a finite energy, while the most discouraging would be the last case, since then the 
cross section would stay exponentially suppressed at all energies. 

The possibility in was advocated [38 - 40] in terms of the so-called "premature 
unitarization" ||0| . The argument is based on considering the interplay in the s-channel 
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Figure 2: Unit arizat ion in the s-channel according to the model[40] of "premature uni- 
tarization". 



unitarity of the processes few — *> many and many — > many. The former processes are 
argued to be still weak (exponentially suppressed) when the processes many — > many 
are at the unitarity limit, which effectively shuts off the further growth of the few — ► 
many cross section. 

A somewhat simplified picture of this behavior is shown in figure 2, where the total 
cross section is represented as imaginary part of a 2 — ► 2 forward scattering amplitude 
through an instanton (J) - antiinstanton (J) configuration. According to the model of 
"premature unitarization" the total amplitude is given by summation over instanton 
- antiinstanton chains iterated in the s-channel, i.e. where all the total energy flows 
through the additional (anti)instantons. Each additional I — I pair brings in the factor 
e -2Si B(E), where the "bond function" B(E) is the multi-boson enhancement of the 
one-instanton-induced cross section observed in || |6], 25]. The summation over the I — I 
chains (figure 2) gives 

B(E)e- s > 



a 



tot 



-Si 



Im tt , (6) 

[l + r 1 {B{E)e-^f\ V ; 

where rj = 0(1) is a rescattering factor. If given by eq.(^), the cross section reaches its 
maximum when B(E)e~ Si = 0(1) and its value at the maximum is of order e~ Si , which 
corresponds to the lower bound of 1/2 for the function F(e). The presented reasoning 
is however oversimplified: it assumes that all the (anti)instantons in the chains have 
same fixed size. Relaxing this assumption leads 0] to a lower bound for F(e), which is 
generally different from 1/2. 

At still higher energies the formula (§) gives a falling cross section. However this 
regime is unphysical: initial particle can shake off energy by emitting one or few hard 



6 



bosons, so that the energy in the collision gets back to the one corresponding to the 
maximum. (Emission of hard bosons suppresses the cross section by a few powers of 
the coupling constant, while the gain in the non-perturbative amplitude is exponential.) 
If indeed the "holy grail" function has a minimum at some energy, this would imply 
that above that energy the process can not be described by semiclassical methods, since 
emission of hard quanta becomes essential. 

It turns out however, that the "premature unitarization" and thus the simple picture 
of the s-channel iteration of instanton-antiinstanton correlations is not mandatory and 
apparently depends on specifics of the theory. The known examples of simplified models, 
where the "holy grail" function is indeed bounded from below by 1/2 are the Quantum 
Mechanical problem with a double well potential [42"|| and the soft contribution to the 



scattering through a bounce f43[, f44f in a (1+1) dimensional model of one real field with 



metastable vacuum |45|. (It has been pointed out[^6|, |47j that in the latter model there 



is also a hard contribution to the bounce-induced scattering, for which the "holy grail" 
function goes to zero at the analog of the sphaleron energy.) Another example, where 
the "holy grail" function is bounded by a value, smaller than 1/2, namely 0.160, is the 
problem of catalysis of false vacuum decay in (3+1) dimensions by collision of two (or 
few) particles f48|. In this problem the semiclassical probability reaches maximum at the 



top of the energy barrier. 

2.2.2 Rubakov - Tinyakov approach. 

The main difficulty in developing a semiclassical approach to the few — > many scattering 
is the presence of hard quanta in the initial state, which state is thus not a semiclas- 



sical one. It has been suggested [(49j p0| to circumvent this difficulty by considering a 
scattering, where a finite small number of particles in the initial state is replaced by 
ni(nitiai) = v/g 2 , where g is the coupling constant in the theory and v is a parameter. 
For a finite v the initial state of this kind can be treated semiclassically, and in the end 
the limit of the probability at v — > 0, or v —>■ const/g 2 is to be considered in order to 
relate to the process few — > many. Within such setting the "holy grail" function depends 
on v: F(e, v) and it is conjectured that its limit at v — > is smooth, which conjecture 
is supported by high-order perturbative calculations around the instanton|51[j. The cen- 
tral point of this approach is that the function F(e, v) is determined from a solution 
to a well-defined boundary value problem[^2j for classical field equations, although in 
essentially complex time. 
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Figure 3: Contour in the complex time plane and boundary conditions for a classical 
solution, describing the scattering of a semiclassical initial state into multiparticle final 
states [52]. 



The classical solution that describes the path of largest probability in a model with 
one real field <fi is evolving along the contour in the complex time plane shown in figure 
3. At Ret — ^ +oo the solution is required to be real, thus its momentum components 
should be of the form 0(k) = b^e' 1 ^ 1 + b*_ k e 1Ukt , while at Ret — > — oo the positive 
frequency part is rescaled by the parameter e e : 0(k) = /k e"'^' + e 9 f^e lbJkt . The 
parameter 9 in the boundary condition and the parameter T of the contour (cf. figure 
3) are Legendre-conjugate of respectively the multiplicity rii and the total energy E of 
initial particles. Namely, if % S is the classical action on the whole contour (thus S is 
defined in the way that it is real in the Euclidean space), one finds |52| 



n, = 2 w , E = 2-. (7) 
Furthermore the "holy grail" function, entering the WKB estimate of the total cross 



section as a tot ~ exp {—g 2 F(e, v)) is given by the Legendre transform of the action: 



2S-ET- nS 



(8) 



Quite naturally, the formulated classical boundary value problem is not easily solv- 
able, and a sufficiently good approximation to the solution is known only in a few 
models[53 - 55]. In particular the model, considered in |5^| , describes one scalar field in 
(1+1) dimensions with the potential 



2 2 2 

m z 2 m z v z 



exp 



(9) 
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Figure 4: The "holy grail" function normalized to F(0, v) = 1 in the (1+1) dimensional 
model [54] with exponential interaction. 



where v and A are dimensionless constants, which both are assumed to be large. The 
parameter 1/v is the small coupling constant of the perturbation theory in the vacuum 
= 0. The negative sign of the interaction term implies that the energy is unbounded 
from below at large 0, thus the vacuum = is metastable, and is separated from the 
decreasing part of the potential by a barrier located at ~ v, provided that A ^> 1. 
Beyond the maximum the potential rapidly goes down, so that the potential essentially 
is a quadratic well with a "cliff" []54f . The metastability of the perturbative vacuum at 
= does not show up in calculations of the scattering amplitudes to any finite order 
of the perturbation theory, and it only arises through a non-perturbative effect: unitary 
"shadow" from the false vacuum decay, which makes this contribution analogous to 
instanton-induced scattering amplitudes in a Yang-Mills theory |45|, |56|. The analog of 
the sphaleron energy is the height of the barrier separating two phases: E$ p = const ■mv 2 . 

At large A the potential ([5]) contains a sharp matching of the quadratic part (free field) 
and a steep exponential "cliff", which enables[54] to solve the boundary value problem 
in the leading order in 1/A and also to clarify the contribution of multi-instanton (multi- 
bounce) configurations. It has been found [54] that the multi-instanton configurations in 
this model are still not important when the one-instanton contribution becomes large. 
As a result the "holy grail" function, as shown in figure 4, reaches zero at finite energy, 
which energy increases when the semiclassical parameter of the initial state multiplicity 
v = riiV 2 decreases. In figure 4 is also shown the behavior corresponding to the periodic 
instanton, which maximizes over rii the rate of tunneling through the barrier in the 



processes n« — > rif at given energy E [1C 
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2.2.3 Prospects for QCD hard processes. 

It has been argued [57 - 60] that a manifestation of instanton- induced scattering in a 
weak coupling regime can be observed in hard processes in QCD. The suggestion is to 
search for final states in hadron collisions, which contain a large number of minijets, each 
with a typical invariant mass /i, such that a s (fi) is sufficiently small, e.g. fi ~ 4GeV, so 
that a s (fi) ~ 0.25. An instanton-induced process should involve production of typically 
rij « 47r/a s (yu) m 50 such jets, which requires energy in a parton - parton collision 
of somewhat higher than rij fi ~ 200 GeV. The prospects of observing the instanton 
induced hard processes in QCD are certainly more phenomenologically attractive, since, 
unlike the electroweak case, the energy range is hopefully within the reach of LHC and 
also the cross section can be of a more encouraging magnitude, even if it is suppressed 
by an exponential factor, like exp(— 2n/a s (fij) ~ 10~ n as suggested by the "premature 
unitarization" models. However the reality of observing these possible non-perturbative 
hard processes in QCD is still under discussion. 



3 Multi-particle production in topologically trivial 
sector 

The growth of the rate of the instanton-induced processes is associated with production 
of multiboson final states until at high multiplicity rif ~ 1/g 2 the final state becomes 
not tractable perturbatively. A similar problem in fact arises [SI, |62| at those high 



multiplicities in processes, which do not require contribution of field configurations with 
non-trivial topology, and thus are allowed in perturbation theory. This is related to the 
well known factorial growth of coefficients in the perturbation theory series|7j. Namely 
in the perturbation theory the total cross section for production of n bosons interacting 
with a weak coupling g is given, modulo the phase space suppression at finite energy, by 

<r n ~ n\ (g 2 ) n . (10) 

At small n, naturally, the cross section is decreasing with multiplicity. However at 
n ~ 1/g 2 the growth of n\ becomes faster than the decrease of (g 2 ) n and the behavior 
(|T0|) would imply that the cross section starts to grow with multiplicity. Therefore the 
question: "If there is enough energy to produce ^> 1/aw W, Z, H bosons, will they 
be actually produced with non-negligible cross section?" does not seem to be entirely 
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paradoxical or idle in view of eg. (|T0D . The difficulty of answering this question is in that 
the lowest order equation ( |T0D becomes inapplicable already at n ~ l/g, since the loop 
corrections to a n are governed by the parameter n 2 g 2 . The latter can be seen from the 
number of rescatterings between the final particles: 0(n 2 ), each having strength g 2 . In 
what follows we will discuss several steps that have been attempted toward answering 
the above question. It also may well be that a solution of the multiboson problem 
without the topological complications will provide an insight into the problem of (B+L) 
violation in high-energy collisions. 

3.1 Multiboson amplitudes in A0 4 theories. 

One of simplest models, where the development of non-perturbative dynamics in multi- 
boson amplitudes can be studied, is that describing one real scalar field with the A 4 
interaction, whose potential is given by 

W) = ^0 2 + ^ 4 - (11) 

If m 2 is positive the field has one vacuum state at (O|0|O) = and the symmetry under 
sign reflection: <ft — > —(f) is unbroken, while at negative m 2 there are two degenerate vacua 
(O|0|O) = ±f with v = |m|/vA, which situation describes the spontaneous symmetry 
breaking (SSB). 

3.1.1 Multiboson amplitudes at zero energy and momentum. 

The simplest problem concerning multiboson amplitudes is, perhaps, that of calculating 
connected n-boson off-shell scattering amplitudes A n , in which all the external particles 
have zero energy and momentum |63|, |64j|. The amplitude A n can be written in terms of 
the connected part of the Euclidean-space correlator: 



J(Ja(x)d d x) exp(-g[0])P0 

fex V (-S[<j>))V<f> ' (12) 

where cr(x) is the deviation of the field <j) from its vacuum mean value cr(x) = <f)(x) — 
(O|0|O), and the integral / a(x)d d x is understood as the p — > limit of the Fourier 
transform / a(x)e ipx d d x. Furthermore the connected part of the correlator ( ]T2|) is con- 
veniently given by the n-th logarithmic derivative of the generating functional Z(j) = 
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/ exp(— S[<f>] + J jcr(x) d d x) V(j) with a constant source j: 

(13) 



Ah 



3=0 



Introduction of a constant source is equivalent to replacing the original potential V(<f>) 
by V(<f>) — j o. Furthermore, for a constant source In Z(j) is related to the energy of the 
vacuum E(j) in the presence of j: lnZ(j) = —VTE(j), where VT is the normalization 
space-time volume. Thus according to eq.(^) the asymptotic at large n behavior of 
the amplitudes A n is related to the position j c of the nearest to j = singularity of 
E(j) in the complex j plane: A n ~ n\j~ n . At the classical level the position of the 
singularity is determined by the value of j, at which two solutions of the equilibrium 



equation dV/dcf) = j coincide. For the potential (|TTD this happens at j c = ±z y'4/27Am 3 , 
which determines the asymptotic behavior |86], EM of the tree-level amplitudes A n : 



l^l -"'(trf)'" 2 - < 14 > 

In a theory with unbroken symmetry the quantum loops modify E(j) according to 
the Coleman- Weinberg potential |>3] thus shifting and modifying the singularity in the 
j plane. However these corrections neither eliminate the singularity nor bring it to 
j = 0. The shift of the position can be absorbed in normalization of A and m, while the 
modification of the type of the singularity only affects sub-leading in n factors, so that 
the leading behavior in eq. ([14]) is not modified by quantum corrections in a theory with 
unbroken symmetry. 

The situation with quantum effects in a theory with SSB is drastically different: 
non-perturbatively the point j = is in fact a branch point of the vacuum energy E(j) 
for either of the vacua. Indeed, if, for definiteness, one choses to consider the amplitudes 
A n in the 'left' vacuum: (O|0|O) = —v with v = |m|/vA, and follows the dependence 
of its energy on j, one finds that this state is stable for real j < and is metastable 
at arbitrarily small positive j. Thus at j > the energy E(j) acquires an imaginary 
part given by the decay rate of the metastable vacuum. In this situation the Taylor 
expansion of E(j) is asymptotic and the coefficients are determined by the decay rate in 



the presence of an infinitesimal positive source term. In this situation the calculation p6 
of the false vacuum decay rate in the thin wall approximation is applicable exactly. Thus 
one can readily find the exact non-perturbative asymptotic behavior of the amplitudes 
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ge n in a theory in d space-time dimensions |p4|| : 

— 

nd \ . I Cd A \ 2 / A 



(2 




m 



4-d 



(15) 



with 



3 d r(d/2) 

22d-l <^-l ^d/2 



2 
d-1 



The factorial behavior in eq. Jig) , if interpreted in terms of loop graphs in perturbation 
theory, corresponds to contribution of graphs with n/(d — 1) loops. 

The considered off-shell amplitudes A n are not physical. However one can draw from 
the described exercise at least two, possibly important, theoretical conclusions about 
multiboson amplitudes: 

• the n! behavior suggested by the tree-level analysis is not necessarily eliminated and 
may even be enhanced in the exact result, and 

• the large n behavior of multiboson amplitudes does not have to be universal and may 
in fact be very sensitive to details of the theory. 



3.1.2 Production of on-shell multiparticle states at and above threshold. 
Tree graphs. 

More physical, than the previously discussed off-shell amplitudes, are the amplitudes 
of processes, where n on-shell bosons are produced by a highly virtual field (f) (1 — > n 
process): a n = (n|0(O)|O). (These e.g. can be related to the reaction e + e~ — > nH.) As 
will be explained, it turns out that one can explicitly find the sum of all tree graphs 
and all one-loop graphs for these amplitudes at any n, provided that the final bosons 
are exactly at rest in the cm. system. Also the summation of two- and higher- loop 
graphs is in principle possible for this kinematical arrangement, however a calculation 
with a finite number of loops is inevitably plagued by the breakdown of the perturbation 
theory at large n. Thus far three methods have been used in calculation of the threshold 
amplitudes of the 1 — > n processes: the Landau WKB method, the recursion equations, 
and the functional technique. 

Landau WKB method |68| is used in Quantum Mechanics for calculating transition 
matrix elements between strongly different levels. (For a field theory derivation of this 
technique see ||69|| .) In the tree graphs for the threshold 1 — > n amplitudes all the 
external and internal lines carry no spatial momentum. Thus the problem is reduced to 
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dynamics of only one mode of the field with spatial momentum p = 0, i.e. to a Quantum 
Mechanical problem. This approach yields the result 170] for the sum of the tree graphs 
at the threshold with accuracy up to terms 0(l/n 2 ) at large n. (Application of the 
Landau WKB technique in the problem of multiboson amplitudes is also discussed in 

Recursion equations for the amplitudes a(n) arise from inspecting the construc- 



tion of Feynman graphs. For the simplest case of tree graphs in A0 theory the algebraic 
form of the equations is 



(V - 1) 



a[n) 



odd ni, 712 



a{ni) a(n 2 ) a(n — n x — n 2 ) 



nil n 2 \ \n — n\ — n 2 ) 



(16) 



where the sum runs over odd n\ and n 2 as well as n is odd, since due to the unbroken sign 
reflection symmetry the parity of the number of particles is conserved. Also the mass 
m in eq.(|l(J) is set to one, since it can be restored in the final result from dimensional 
counting. The solution to the equations (|T6| ) reads as[|64[ 



a{n) = n\ (A/8m 2 )" 



(17) 



which can be found by applying the regular method of generating functions||73[|. For the 
theory with SSB the recursion equations are modified by the presence of cubic vertices. 
The result for the amplitudes in the theory with SSB is|73| 



a\n) 



-n\ {2v 



,1-n 



18) 



The recursion method can be extended to other theories |74| as well as to loop 
graphs|75|, |76[ and to an analysis of higher loops|77|]. However a more convenient method 



for further analysis is the one suggested by Brown |78| and is based on a functional tech- 
nique. Before proceeding to discussing this method and its further applications we report 
on estimates of the tree amplitudes above the threshold and thus of the total probability 
of the processes 1 —>■ n at a high energy E. 



3.1.3 Lower bound for cross section at the tree level. 



The tree graphs for the processes 1 — > n in a A0 theory all have the same sign EL], |62 



The decrease of the amplitude above the threshold is thus determined by the increasing 
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virtuality of the propagators in those graphs, which depends on the kinematics of the 
final state. One can thus find a lower bound on the tree amplitudes above the threshold 
in a restricted part of the final-particle phase space[79 - 81], which gives a lower bound 
on the total probability of the process. In particular, if the kinematical restriction is 
chosen [J7S|, [SD| as the condition that the cm. energy of each individual particle in the 
final state does not exceed lu, then in this region of the phase space the tree amplitude 
A(l — > n) is larger than the threshold amplitude a(n) in which the physical mass M of 
the scalar boson is replaced by |u; in the theory without SSB (in which case M = m) 
and by |o> in the theory with SSB (where M = \/2 \m\). The cut off energy uj is then 
optimized for each value of the total energy E and multiplicity n in order to find the 
largest lower bound on the total probability 

a n = J \A(1 -> n)\ 2 dr n , (19) 

which is given by the integral over the n particle phase space r n . As a result the lower 



bound on a n is found |8(| in the scaling form 

"47T 2 C 



a n > exp 



A 



(20) 



where v = nM/E is the ratio of the multiplicity n to its maximal possible value E/M, 
e = E / Eq with Eq being an analog of the 'sphaleron' energy: E = 4vr 2 cM/A, and 
the constant c in these formulas is c = 9 (c = 8/3) for a theory without (with) SSB. 
The calculated |8(| behavior of the function f(e, v) is shown in figure 5, which thus 



illustrates and quantifies the interplay between the n\ and the power of small coupling 
constant, discussed in connection with eq.(|l(]). The function f(e, v) displays a normal 
perturbative maximum at zero multiplicity. However, as energy grows, and production 
of high multiplicity states becomes unsuppressed kinematically, this function develops 
a second maximum, which at larger energies eventually crosses zero with an apparent 
violation of unitarity. 

It is interesting to notice that the kinematical suppression of multiparticle final states 
is quite essential and shifts the energy, at which the tree graphs violate unitarity signifi- 
cantly higher than one would guess from a simple estimate E cr n ~ 47r Mj A. If applied to 
a multi-Higgs production in the Standard Model the lower bound (^) breaks unitarity 
at E mt w 15.5 (32n 2 M H /\) w 1000 GeV(200GeV/M H ). 

It should be also mentioned that it is quite likely that the scaling behavior (|20|) at a 
given large multiplicity n also holds for the actual cross section, which point is recently 
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0.2 0.4 0.6 0.8 1 

v 

Figure 5: The function /(e, v) vs. v at several characteristic values of energy: low energy, 
no secondary maximum (e = 1), the secondary maximum just developed (e = 10), the 
secondary maximum becomes global and is just above the unitarity limit (e = 15.5) [80]. 



strongly emphasized in [82]. The function f(e,is) can thus be called differential in n 
"holy grail" function. 



3.2 Generating field technique. 
3.2.1 Tree level. 

A more convenient and more conceptually transparent technique for dealing with tree- 
level threshold multiboson amplitudes was suggested by Brown [fT8| and was later ex- 
tended to calculation of one-loop |83|, 



and higher quantum effects [g5|, |86[ in these 
amplitudes. The technique is based on the standard reduction formula representation 
of the amplitude through the response of the system to an external source p(x), which 
enters the term p<j) added to the Lagrangian. 



n\cj)(x)\Q) 
6 



Pi) 



5p(x a ) 



n . 

TT lim / d^Xa e WaXa * 
(0ouMx)\0 in )0\ p=0 , 



m 



(21) 



the tree-level amplitude being generated by the response in the classical approximation, 
i.e. by the classical solution (po(x) of the field equations in the presence of the source. 
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For all the spatial momenta of the final particles equal to zero it is sufficient to 
consider the response to a spatially uniform time-dependent source p(t) = p (u) e luJt 
and take the on-mass-shell limit in eq.(|2~T|) by tending uj to m. The spatial integrals 
in eq.(EHp then give the usual factors with the normalization spatial volume, which as 
usual is set to one, while the time dependence on one common frequency uj implies that 
the propagator factors and the functional derivatives enter in the combination 

v Fa! 5p{x a ) v ! 5p{t) 5z(t) v ; 

where 

zit) = 23 

coincides with the response of the field to the external source in the limit of absence 
of the interaction, i.e. of A = 0. For a finite amplitude po of the source the response 



z(t) is singular in the limit to — > m. The crucial observation of Brown|]78[ is that, since 
according to eq. Q2"2"D we need the dependence of the response of the interacting field (f> 
only in terms of z(t), one can take the limit po{uj) — > simultaneously with u — > m in 
such a way that z{t) is finite: z{t) — > Zoe tmt . 

Furthermore, to find the classical solution (j>o(x) in this limit one does not have to 
go through this limiting procedure, but rather consider directly the on-shell limit with 
vanishing source. The field equation with zero source in A0 4 theory without SSB is 

d 2 (j) + m 2 (f) + A0 3 = . (24) 



For the purpose of calculating the matrix element in eq.(|2"T|) at the threshold one looks 
for a solution of this equation which depends only on time and contains only the positive 
frequency part with all harmonics being multiples of e tmt , which condition is equivalent 
to requiring that (f>(t) — > as Imt —> +oo. The solution satisfying these conditions reads 
as[[F|] 

^ = i-ww)' (25) 

According to equations (|2~2|) and (^l|) the n-th derivative of this solution with respect 
to z gives the matrix element (n|0(O)|O) at the threshold in the tree approximation: 

2fc+l 



(2fc+l|0(O)|O) o 
which reproduces the result in eq. (|17D . 



dz 



(26) 

z=0 
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Figure 6: The tadpole graph in the background field O for calculating one-loop cor- 
rection to 0o- The double lines and the heavy dot represent respectively the Green's 
function and the vertex in the background field. 



For the case of theory with SSB the solution reads as 



which reproduces the tree amplitudes in eq. fll~8|) . In this case z(t) = e lMt , where M = 
v2 \m\ is the mass of physical scalar boson. 

3.2.2 One-loop level. 

To advance the calculation to the one-loop level one has to calculate the first quantum 



correction (j)\{t) to the classical background field 0o- This amounts [|S3|1 to evaluating the 
tadpole graph of figure 6, where both the Green's function and the vertex are calculated 
in the external background field 0o- The green's function G(x; x') satisfies the equation 

(d 2 + m 2 + 3A 0o (*) 2 ) G{x; x') = -i 5{x - x') , (28) 

in which the differential operator in the Minkowski time contains explicitly complex 
field 0o (cf. eq. (|25|) or eq.(p7|)). A straightforward rotation to the Euclidean time, 
i t — > r, is problematic, since the background field then develops a pole at a real r. The 
acceptable solution is achieved by simultaneously rotating and shifting the time axis in 
eq.(|28D in such a way that —Xz(t) 2 /8m 2 — > exp(2mr) for the theory without SSB, and 
—z(t)/2v — > exp(Mr) for the theory with SSB. In terms of thus defined r the background 
field has the form 0o(t) = i ^2/\m/ cosh(mr) (no SSB) and 0o(t) = v tanh(Mr/2) 
(with SSB). In both cases the term 4>o(t) 2 in equation (^8|) is real and non-singular. After 
applying the standard decomposition of the Green's function over the conserved in the 
background 4>o(t) spatial momentum k: 

G(r, x; t', x') = / G w (t, t') , (29) 
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one arrives for the case of no SSB at the well-known in Quantum Mechanics equation 

/ + - 2 " jJL^) = S(r ~ r>) (30) 



dr 2 (cosh t 



with uj = Vk 2 + 1 and the mass m set to one. (For the theory with SSB one gets the 
same equation with a rescaled o;||84||.) Thus the problem of finding the first quantum 
correction (f>x(t) to the background field is completely solvable on the r axis, and the 
solution can then be extended to the whole complex plane of t by analytical continuation. 
For the theory with no SSB the result |83|] for the amplitudes a{n) at the one- loop level 
reads as 

a(n)o+i = a(n) (l-(n- l)(n - 3) ^ Fj , (31) 

where 

F = 4L 2 -±4-i.) . (32) 
The analog of this result for the SSB case is[]84| 



a(n) +i = o(n) ^1 +n(n - 1) • (33) 

3.3 Nullification of threshold amplitudes. 

The equations (|31|) and ( p3[ ) display a remarkable feature: in spite of the presence 
of an intermediate state with two bosons in one-loop graphs, their contribution to the 
amplitudes in the case of SSB is real, while the factor F in eq.(^) is an easily recognizable 
threshold factor for the 2 — > 4 process with no indication of presence of other thresholds. 
Using the unitarity relation for the imaginary part of the loop graphs, one immediately 
concludes that this can only be if the tree amplitudes of the on-shell processes 2 — > n 



are all zero at the threshold for n > 4 in the theory without SSB|83], |87| and for n > 2 
in the theory with SSBQ . 



This can be traced to the special properties of the reflectionless potential — 6/(cosh r) 2 
in equation ( p0|) and generalized |88|, [T6| to other theories, where the problem of the 



2 — ► n scattering is reduced to finding the Green's function in the reflectionless potential 
—N(N + l)/(coshr) 2 with integer N. The known additional cases are the following: 

• Linear a model (N = 1): the tree-level amplitudes of the scattering it it — > na are 
all zero at the corresponding thresholds for n > 1 . 
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Fermions with Higgs-generated mass: if 2m//mfl = N (integer), then all tree- level 
amplitudes of / / — > n H are zero at threshold for all n > N. 

Vector bosons with Higgs-generated mass: if Arriy/mfj = N(N + 1), then all tree- 
level amplitudes with transversal vector bosons of Vr Vr n H are zero at thresh- 
old for all n > N. For longitudinal vectors the same amplitudes of Vj, Vi — > n H 



are zero for n = N [8J| and for n > N + 1 [76|. 

All these cases, except the one with longitudinal vectors, stem[^] from the generic 
interaction of two fields and x °f the form | 4> 2 x 2 and the self-interaction of the field 
4> as described by the potential (|H]). Then if the ratio of the coupling constants satisfies 
the relation 2£/A = iV(iV-|- 1) with N integer, the tree-level threshold amplitudes of the 
processes 2% — > ncf) are all zero for n > N in a theory with SSB and for n > 2N in a 
theory with unbroken symmetry. This behavior somewhat resembles the nullification of 
inelastic amplitudes in the Sine-Gordon theory, where it is a consequence of a symmetry 
and is a deep property of the theory. In the theories considered here this is a much 



weaker property, which holds only at threshold and, generally, only at the tree level [90 



However the nullification in this case can be a consequence of a hidden symmetry, which 
holds at the classical level and/or has a more restricted scope. Thus far such symmetry 



has been revealed [91] only for the case of N = 1, where it can be traced to the symmetry 
of a system of two anharmonic oscillators, described by the potential 

V(x,y) = ^ + ^ + ^ + yr. (34) 

If the frequencies u)\ and io 2 were equal, the model would have an 0(2) symmetry, 
corresponding to conservation of the angular momentum Q = xy — xy. However even 



for lui ^ co>2 the symmetry persists ]^T|] corresponding to conservation of the invariant 

f q 2 + K 2 - c 2 )^ 2 + f y 2 + \tf + f *V). 

It should be also noted that if the ratio 2£/ A does not satisfy the above mentioned 
condition, the threshold amplitudes of the processes 2% — » ncj) display|^| a 'normal' 
factorial growth with n. 



3.4 Non-perturbative analysis. 

The n 2 \ behavior of the loop corrections in the equations (|3lD and (|33|) convinces us 
that the perturbation theory is of little help in finding the amplitudes at large n and 
a true non-perturbative analysis is required. It turns out that to a certain extent such 
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analysis can be performed for the A0 4 theory with SSB. In terms of the variable r the 
problem of calculating the threshold amplitudes a(n) = (n|0(O)|O) reduces [55] to a well 
defined Euclidean-space problem of calculating the quantum average $(r) of the field 

$(T) = Je-W0 (35) 

with the kink boundary conditions, i.e. — > ±u as r — > ±00. The average field then 
expands at r — > —00 in the series 

CO 

Hr) = J2c n e nMT (36) 

n=0 

and the threshold amplitudes are given by a(n) = n\c n /ci, where the coefficient c\ 
describes the one-particle state normalization: c\ = (1|0(O)|O). 

Due to the fact that the classical kink solution provides the absolute minimum for 
the action under specified boundary conditions, the path integrals in eq.(p5[) are well 
defined (no negative modes) and thus the average field $(r) is real in any finite order 



of the perturbation theory. Thus the coefficients c n of the expansion fl3~6|) are real too. 
Therefore the amplitudes a(n) are real to any finite order in A. As we have seen at the 
one loop level this implies the nullification of the tree amplitudes of 2 — > n for n > 2. In 
higher loops this implies a relation between the amplitudes of the processes k — > n with 
different k. The only exception is the particular case of n = 3, for which the imaginary 
part of the a(3) can be contributed only by the two-boson intermediate state. Since the 
imaginary part is vanishing, one concludes ||| that the 2 — > 3 amplitude is vanishing at 
the threshold in all orders in A. 

The function $(r) given by the expansion (|36|) is manifestly periodic: $(r+2z7r/m) = 
$(r). Using this property and the boundary conditions at r — ► ±00 one finds||| that 
the exact function $(t) necessarily has a singularity at a finite r. Thus the expansion 
(p6|) has a finite radius of convergence, and thus the exact threshold amplitudes a(n) 
grow at least as fast as n\. In other words, the quantum effects do not eliminate the 
factorial growth of a(n). 

As is indicated by the n 2 X parameter of the perturbation theory for the coefficients c n , 
the saddle point (SP) for the action S [(/)], given by the x-independent 'domain wall' with 
the kink profile is not the correct SP for calculating the coefficients c n at large n. It has 
been argued [j86f that the correct SP configuration is given by a spatially inhomogeneous 
field configurations in which the 'domain wall' is deflected towards negative r by a 
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maximal amount ho. Then at a large negative r the coefficients c n are given by 

c n ~ e nMho e'^ A - Ao) , (37) 

where A is the area of the domain wall with deflection, A is the same for the undistorted 
flat wall, and \x = |m| 3 /3A is the surface tension of the wall. 

Finding the extremum of the 'effective action' in the exponent in eq.(|3"7|) exactly 
corresponds to finding the equilibrium configuration of a (d — 1) dimensional membrane 
with a force equal to nM applied at the point of maximum deflection. In general this 
problem has no real solution (the film gets punched). However a solution exists, where 
a part of the trajectory of the domain wall resides in the Euclidean space, and a part 



is in the Minkowski space |pq| . The Minkowski-space part of the trajectory corresponds 
to evolution of a bubble made of a domain wall and having energy E = n M . The 
amplitudes a(n) are then found as a sum of resonant contributions of the quantized 
levels of the bubble: 

a(n) ~n\ i _ exp [f(d) E (E , (38) 

where f(d) is a positive coefficient depending on the space-time dimension, and 1(E) = 
§ pdr is the action of the bubble over one period of oscillation. Clearly the amplitudes 
a(n) in eq . (|38D have poles at the values of E satisfying the Bohr- Sommerf eld condition 
1(E) = 2irN. 

The result in eq.(^) can be interpreted as that the growth of a(n) is due to a strong 
coupling of the bubble states to the multi-boson states with all particles being exactly at 
rest. However, it is known from the numerical studies of mid-70s[93 - 95] that the lifetime 
of the bubbles is of order one in units of their period. Thus one should conclude |86|] that 
there arises a non-perturbative form factor, which cuts off the integral over the phase 
space of the final bosons and makes the total probability of a moderate value, inspite of 
the extremely large value of the coupling to exactly static bosons. The total probability 
of the process 1 — > n in this picture is given by the probability of creating a bubble with 
energy E by a virtual field 0. This can be estimated ||86|| by the Landau WKB method 
and is found to be exponentially small: 

a(l -> B(E)) ~ exp [-2f(d)E(E//j)^] , (39) 

where f(d) is the same as in eq . (|38D . Thus one concludes that in the theory with SSB the 
total cross section of non-perturbative multiparticle production is extremely likely to be 
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exponentially small at high energy. However because of the usage of special properties 
of the theory with SSB it is not clear, whether this conclusion can be generalized to 
other theories, in particular, to the multi-boson production in the Standard Model. 

4 Conclusions. Problems 

The problem of multi-particle processes in theories with weak interaction is one of most 
challenging in the quantum field theory. In solving this problem we are most likely to 
find new methods of non-perturbative analysis of the field dynamics. As it stands now, 
there are mostly problems facing us, some of which are: 

• It is not clear, to what extent the exponential suppression of the (B+L) violation 
in particle collisions is lifted at high energy: by a factor 1/2 in the exponent, by a 
different factor, or completely. All these types of behavior are observed in simplified 
models. 

• The n\ behavior of the amplitudes for production of n bosons survives the quan- 
tum effects, at least in some models. However this does not necessarily imply a 
catastrophic growth of the cross section. 

• Peculiar zeros are observed in threshold amplitudes of multi-boson production. 
However it is not clear, whether they signal some deep properties or this is a mere 
coincidence. 

• The classical field configurations give rise to multiparticle amplitudes. However 
their role in high-energy collisions is yet to be understood. 

This work and the author's participation in the conference are supported, in part, by 
the DOE grant DE-AC02-83ER40105. 
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Discussion 

A. Kataev, CERN and INR, Moscow: 

Do the theoretical results discussed in your talk have any interesting phenomenological 
implications? Is it possible to study these non-perturbative effects experimentally? 

M. Voloshin: 

Possible effects in hard processes in QCD with production of many minijets might be 
observable at LHC energies. An experimental study of possible non-perturbative pro- 
cesses in the electroweak interactions may require an energy of about 1000 TeV. 
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V. Kuvshinov, Minsk: 

It seems we have here new mechanisms for multiparticle production. For example, it 
can give contributions to the intermittency phenomenon. Is multiplicity important for 
B + L violation? Or is only n\ important? 

M. Voloshin: 

The multiplicity is important through the n\, or, possibly, a stronger factor. 
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